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; ABSTRACT 

'. This paper describes Herschel observations of the nearby (8.5pc) G5V multi-exoplanet 

■ host star 61 Vir at 70, 100, 160, 250, 350 and 500 j^m carried out as part of the 

CnI ■ DEBRIS survey. These observations reveal emission that is significantly extended out 

to a distance of > ISarcsec with a morphology that can be fitted by a nearly edge-on 
(77° inclination) radially broad (from 30AU out to at least lOOAU) debris disk of 
, fractional luminosity 2.7 x 10~^, with two additional (presumably unrelated) sources 

. nearby that become more prominent at longer wavelengths. Chance alignment with 

a background object seen at 1.4GHz provides potential for confusion, however the 
star's 1.4arcsec/year proper motion allows archival Spitzer 70/im images to confirm 
that what we are interpreting as disk emission really is circumstellar. Although the 
exact shape of the disk's inner edge is not well constrained, the region inside 30AU 
must be significantly depleted in planetesimals. This is readily explained if there are 
additional planets outside those already known (i.e., in the 0.5-30AU region), but 
is also consistent with collisional erosion. We also find tentative evidence that the 
presence of detectable debris around nearby stars correlates with the presence of the 
lowest mass planets that are detectable in current radial velocity surveys. Out of an 
unbiased sample of the nearest 60 G stars, 11 are known to have planets, of which 6 
(including 61 Vir) have planets that are all less massive than Saturn, and 4 of these 
have evidence for debris. The debris toward one of these planet-hosts (HD20794) is 
reported here for the first time. This fraction (4/6) is higher than that expected for 
nearby field stars (15%), and implies that systems that form low-mass planets are 
also able to retain bright debris disks. We suggest that this correlation could arise 
because such planetary systems are dynamically stable and include regions that are 
populated with planetesimals in the formation process where the planetesimals can 
remain unperturbed over Gyr timescales. 

Key words: circumstellar matter - planets and satellites: formation - planet-disc 
interactions. 
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I INTRODUCTION 

Main sequence stars, like the Sun, are often found to be 
orbited by circumstellar material that can be categorised 
into two groups, planets and debris, the latter of which is 
comprised of asteroids, comets and the dust derived from 
them (e.g., Wyatt 2008; Krivov 2010). Although there are 

II examples of nearby stars that are known to have both 
planets and debris (e.g., Moro-Martm et al. 2010), there is 
as yet no evidence for any correlation between the two types 
of material (Greaves et al. 2004; Kospal et al. 2009), and the 
properties of the debris disks around stars that have planets 
are not found to be significantly different to those of stars 
without known planets (Bryden et al. 2009). This is usually 
explained as a consequence of the spatial separation between 
the planets (which are typically found within a few AU) and 
debris (which is typically found at 10s of AU). 

Despite this spatial separation, it is still expected that 
outer debris can be affected by the gravitational perturba- 
tions of close-in planets (Mustill & Wyatt 2009), and that 
such debris can have a significant dynamical effect on the 
evolution of interior planets, e.g., through promoting migra- 
tion (Kirsh et al. 2009) or by triggering instabilities (Tsi- 
ganis et al. 2005; Levison et al. 2011). The delivery of de- 
bris to the inner planets may also be an important source 
of volatiles to the planets (e.g., Horner & Jones 2010). Fur- 
thermore it is reasonable to expect that the initial conditions 
in protoplanetary disks that favour the formation of certain 
types of inner planetary system might also result in (or ex- 
clude) specific types of debris disk (e.g., Wyatt, Clarke & 
Greaves 2007; Raymond et al. 2012). Thus the search for 
any correlation between the two phenomena continues (e.g., 
Maldonado et al. 2012), so that light may be shed on the 
interaction between planets and debris, as well as on the 
formation mechanism for the two components. 

This paper describes observations carried out as part 
of the key programme DEBRIS (Disc Emission via a Bias- 
free Reconnaissance in the Infrared/Submillimetre) on the 
Herschel Space Observatory^ (Pilbratt et al. 2010). DEBRIS 
is an unbiased survey searching for evidence of circumstel- 
lar dust at 100 and 160 /im toward the nearest ~ 80 stars 
of each spectral type A,F,G,K,M (see Phillips et al. 2010 
for a description of the sample). The first results have al- 
ready shown that Herschel observations have the potential 
to detect disks down to much fainter levels than previous ob- 
servations, and moreover have the resolution to resolve the 
disks at far-IR wavelengths (Matthews et al. 2010; Churcher 
et al. 2011; Kennedy et al. 2012). 

Several of the stars in the sample are known planet 
hosts; others may be found to host planets in the future. As 
such, the survey is ideally suited to determining whether any 
correlation exists between planets and debris. Here we focus 
on Herschel observations of the star 61 Vir (HD 115617), 
which is a main sequence G5 star. At 8.55 di 0.02 pc (van 
Leeuwen 2007) 61 Vir is the 8th nearest main sequence G 
star to the Sun. This star exemplifies the merits of an un- 
biased survey, since the star was relatively unremarkable 
at the time the survey was conceived, but has since been 

^ Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and 
with important participation from NASA. 



found by radial velocity surveys to host 3 planets at 0.05, 
0.218, 0.478AU with minimum masses (i.e., Mpisini) of 5.1, 
18.2, 22.9 M© respectively (Vogt et al. 2010), and has sub- 
sequently been the subject of several studies on the secu- 
lar dynamics of its planetary system (Batygin & Laughlin 
2011; Greenberg & van Laerhoven 2012) and on its forma- 
tion mechanism (Hansen & Murray 2012). 

Thus 61 Vir is one of the first of a growing number 
of systems around which only low-mass planets are known, 
which we define here as sub-Saturn mass. Such low-mass 
planets have only been discovered recently, often in multi- 
ple planet systems, either by high precision radial velocity 
measurements (Lovis et al. 2006; Mayor et al. 2011), or by 
transiting studies with Kepler (Borucki et al. 2011; Lissauer 
et al. 2011). 

Previous Spitzer observations showed that 61 Vir hosts 
a debris disk with a fractional luminosity of around 2 x 
10"^ (Bryden et al. 2006), but the location of the emission 
was uncertain, with estimates ranging from 8.3AU (Trilling 
et al. 2008) and 4-25 AU (Lawler et al. 2009) to 96-195AU 
(Tanner et al. 2009). Our observations confirm the existence 
of a bright debris disk, and moreover the resolution of the 
images permit the disk structure to be ascertained. This 
allows us to consider the relationship between the disk and 
the planetary system, and the implications for the system's 
formation and evolution. Given that there are other systems 
with debris in low-mass planetary systems (e.g., Beichman et 
al. 2011), we also consider whether this discovery heralds an 
emerging correlation between the presence of debris and low- 
mass planets (e.g., as predicted by Raymond et al. 2011). 

The layout of the paper is as follows. The observations 
are described in §2, including both the Herschel observa- 
tions and archival Spitzer and radio observations that are 
needed to disentangle circumstellar from extragalactic emis- 
sion. Modelling of the Herschel observations is presented in 
§3 to derive the radial distribution of dust. We find that the 
known planets probably do not play a primary role in stir- 
ring the debris, and discuss the implications of the debris for 
the evolution of the planetary system in §4. Finally a statis- 
tical analysis of an unbiased sample of nearby stars is given 
in §5 to consider the possibility of a disk-planet correlation. 
Conclusions are presented in §6. 

2 OBSERVATIONS 
2.1 Herschel 

2.1.1 PACS 

Herschel DEBRIS observations of 61 Vir were taken on 2010 
August 10. These consist of simultaneous 100 and 160 /im 
images taken with PACS (Photodetector and Array Camera 
& Spectrometer, Poglitsch et al. 2010), using the small scan- 
map mode. Scan-map observations have eight repeats in two 
scan directions that differ by 40° , at a scan rate of 20 arcsec 
s~^. Eight 3-arcmin scan legs were performed per map with a 
2-arcsec separation between legs. The total observing time 
was 890s. An identical observation was obtained on 2011 
July 3, this time with simultaneous 70 and 160 /im images. 
The Herschel observations are summarised in Table 1. 

The data were reduced using a near-standard pipeline 
with the Herschel Interactive Processing Environment 
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Table 1. Herschel observations of 61 Vir. 20 
Obsid Date Instrument Duration 



1342202551 


2010 August 10 


PACS 100/160 


445s 


1342202552 


2010 August 10 


PACS 100/160 


445s 


1342212412 


2011 January 9 


SPIRE 250/350/500 


721s 


1342223608 


2011 July 3 


PACS 70/160 


445s 


1342223609 


2011 July 3 


PACS 70/160 


445s 



(HIPE, version 7.0 build 1931, Ott 2010). The data were 
pre-filtered to remove low-frequency (1/f) noise using a box- 
car filter with a width of 66 arcsec at 70 and 100 /xm and 
102 arcsec at 160 /xm. Different filter widths were tested 
to confirm that this did not significantly affect conclusions 
about disk fiuxes or morphology. The point spread function 
(PSF) of the PACS instrument includes power on large scales 
(e.g., 10% beyond 1 arcmin) that is removed by the filter- 
ing process. For point sources this can be readily accounted 
for in the calibration by comparison of the aperture fiuxes 
of bright stars with known fiuxes. For extended sources like 
61 Vir, the calibration is slightly more complicated, as de- 
scribed below. 

The pixel scale in the maps was set to 1 arcsec at 70 and 
100 /xm and 2 arcsec at 160 /xm. Since this is smaller than the 
native pixel scale, noise in neighbouring pixels is correlated, 
and to derive the uncertainty on any given measurement 
(such as fiux within an aperture) the noise per pixel in the 
resulting map must be interpreted appropriately (Fruchter 
& Hook 2002; Kennedy et al. 2012). Because the DEBRIS 
survey comprises many stars, we can confirm that this cor- 
rection has been applied correctly, because a histogram of 
(observed-predicted) / uncertainty for all observations is con- 
sistent with a Gaussian of unity dispersion as expected. If 
the histogram had been wider, for example, it would indicate 
that the uncertainties had been underestimated. 

Fig. 1 shows the 70, 100 and 160 /xm images cropped 
to the region ±40 arcsec from the star. The peak emission 
at 70 and 100 /xm was found respectively to be offset (in 
different directions) by 1.8 and 1.6 arcsec from the nomi- 
nal location of the star (13h 18m 23.52s, -18d 18m 51.6s). 
Since this is consistent with the expected pointing uncer- 
tainty of Herschel of 2 arcsec (Icr), a value that in turn is 
consistent with the pointing offsets seen for stars detected in 
the DEBRIS sample, we applied a pointing correction on the 
assumption that the star was at this peak location. As the 
two 160 /xm images were taken simultaneously with those at 
70 and 100 /xm, this re-registering process also determines 
the centre of the 160 /xm maps that were then coadded to 
give the presented image. 

The images at all three PACS wavelengths show a very 
similar morphology, of emission that is peaked on the star, 
but that is extended symmetrically about the star along a 
position angle of ~ 65°. This suggests that the star is sur- 
rounded by a disk that is inclined to our line of sight. After 
subtracting a PSF scaled to the stellar photospheric emis- 
sion (see Fig. 2 for the resulting image at 70 /xm), elliptical 
Gaussian fits to the 3 images find for the FWHM of the ma- 
jor and minor axes and for the position angle: 15.9 ± 0.36", 
7.7±0.2", 64±1.3° at 70 /xm, 17.4±0.7", 8.7±0.3", 65=bl.6° 
at 100 /xm, and 24.8 di 2", 14.8 ± 1", 75 ± 4.5° at 160 /xm. 
Given that the FWHM of the PSF is 5.6, 6.7 and 11.4 arcsec 
at these wavelengths, the emission is extended along both 
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Figure 2. Star-subtracted 70 /xm PACS image of 61 Vir. Con- 
tours are the same as in Fig. 1. 



the major and the minor axes. Quadratic subtraction of the 
FWHM of the PSF suggests a disk inclination of around 
20 ± 2° from edge-on at 70 and 100 /xm and 26 ± 7° at 
160 /xm (see §4.2 for a discussion of the implications for the 
planet masses) ; the fact that the long wavelength emission is 
more extended (after this quadratic subtraction) than that 
at 70 and 100 /xm suggests that the disk may be radially 
broad (since longer wavelengths tend to probe cooler mate- 
rial). The uncertainty in the level of photospheric emission 
to be subtracted from the images has a negligible effect on 
the inferred disk structure; e.g., the resulting uncertainties 
in the parameters of the elliptical Gaussian fit are smaller 
than those quoted above. 

The total fiuxes were calculated from the images using 
aperture photometry. A circular 20 arcsec (170AU) radius 
aperture was employed since, at least for 70 and 100 /xm, it 
was found that there was no fiux present on larger scales. 
The fiuxes derived in this way need to be revised because a 
significant fraction of the fiux has been filtered out in the 
reconstruction of the map. This correction factor has been 
derived for point sources in the DEBRIS survey to be 16, 
19, 21% at 70, 100 and 160 /xm respectively, based on how 
the inferred fiuxes compare with those predicted from the 
photosphere models for those stars for which purely pho- 
tospheric emission is confidently detected, and we consider 
the uncertainty in our correction factor to be ^ 5% (i.e., 
the correction factor at 70 /xm is 16 =b 5%). These values 
are consistent with the aperture correction for a 20 arcsec 
aperture quoted by the PACS team^ of 17, 20 and 26%. Al- 
though this correction factor was derived for point sources, 
it should also apply to resolved sources, as long as the aper- 
ture size is larger than the extended source, a condition that 
is tested for the 61 Vir maps in §3. The derived fiuxes are 
198 ± 3 mJy at 70 /xm, 152 ± 6 mJy at 100 /xm, and 133 ± 6 

^ http:/ /herschel. esac.esa.int/twiki/bin/view/Public 



4 M. C. Wyatt et al 




40 



20 



40 20 -20 

(a - a.2oii)*cos(8) (") 



40 



20 



I 



-20 I 



-40 



-20 



-40 



40 



20 




-20 



-40 



40 



20 -20 

(a - a.2oio)*cos(8) (") 



-40 



40 



40 



20 



20 -20 

(a - a*2oio)*cos(6) (") 




-20 



-40 



40 



20 



-20 



-40 



20 -20 

(a - a*2oio)*cos(8) (") 



-40 




40 



20 -20 -40 

(a - a.2oio)*cos(8) (") 




20 -20 

(a - a.2oio)*cos(8) (") 



-40 



Figure 1. DEBRIS PACS and SPIRE images of 61 Vir, all centred on the star at the epoch of the observation, the location of which 
is shown with an asterisk: (top left) 70 /im, (top middle) 100 /im, (top right) 160 /im, (bottom left) 250 /im, (bottom middle) 

350 /im, (bottom right) 500 /im. Contours shown in black are those of the presented images at levels of 4, 13, 22, 31, 40, 49, 58 times 
(770 = 0.028 mJy/pixel at 70 /im, 4, 9, 14, 19, 24, 29 times crioo = 0.029 mJy/pixel at 100 /im, 4, 8, 12, 16, 20, 24 times crieo = 0.052 
mJy/pixel at 160 /im, 3, 4, 5, 6, 7, 8 times (7250 = 4.8 mJy/beam at 250 /im, 2, 3, 4, 5, 6 times (7350 = 6.2 mJy/beam at 350 /im, and 2, 
2.5, 3, 3.5 times (7500 = 6.4 mJy/beam at 500 /im, where ax is the pixel-to-pixel uncertainty in the images (which is measured from the 
images and so includes confusion noise), and the colour scale goes from —2(7\ to the maximum brightness. 



mJy at 160 /im (noting that these uncertainties need to be 
increased due to the systematic effects outlined in the next 
paragraph) . 

The uncertainties in the previous paragraph have been 
derived from the distribution of fluxes obtained by plac- 
ing apertures of the same size at different locations in the 
map. Given the uneven coverage of the map resulting from 
the observing strategy, the centres of these apertures were 
constrained to lie in regions with uncertainties within 10% 
of that of the region around 61 Vir. Additional uncertain- 
ties of 1.4, 1.6 and 3.5% for repeatability in the different 
wavebands^, and the 5% calibration uncertainty mentioned 
above, also need to be applied. A quadratic combination of 
uncertainties results in ±11, ±10, ±17 mJy at 70, 100 and 
160 /im respectively. 

These fluxes can be compared with those determined 
by MIPS. The 70 /im MIPS bandpass is very similar to 
that of PACS, and the MIPS flux has been quoted as 



^ http:/ /herschel.esac.esa.int/twiki/bin/view/Public 



185.6 ± 16.6 mJy (Bryden et al. 2006) and 195 ± 8 mJy 
(Trilling et al. 2008), but since the source is resolved the 
flux is expected to be higher than early estimates, with 
224 ± 8 mJy quoted most recently (Lawler et al. 2009). In 
§2.2.2 we reanalyse the MIPS data to give 237 ± 10 mJy 
in a 36 arcsec radius aperture. Although the PACS flux is 
lower than that measured by MIPS, the values are consistent 
within plausible statistical and systematic errors (e.g., cross- 
calibration between the instruments). The 160 /im MIPS 
flux is given as between 89 and 141 ±20 mJy (Tanner et 
al. 2009), where different values are obtained depending on 
how the observations are corrected for a spectral leak. These 
fluxes are very close to that measured by PACS, and given 
the issues with the leak, the MIPS flux is not used in the 
subsequent analysis. 

There are two additional lower level features in the 
PACS images. There is clearly a source offset ~ 27 arcsec 
to the west of the star. This is not detected at 70 /im (a 
3cr upper limit of 6m Jy), it is present above the 3cr level 
at 100 /im (flux of 10 ± 2 mJy), and is more confidently 
detected at 160 /im (flux of 19 ± 2 mJy). Given the low 
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signal-to-noise, there is no evidence that this emission is ex- 
tended, and we attribute it to a background source. There 
is also a protrusion to the north of the star which occurs in 
the 4(7 contour at 70 and 100 //m, and in the 12a contour 
at 160 /xm. For the shorter wavelengths this level of pertur- 
bation is consistent with that expected from noise superim- 
posed on a symmetrical disk, but the deviation at 160 /xm 
cannot be explained in this way. We considered whether this 
is caused by the shape of the PACS beam, which at such low 
levels is tri-lobed (noting that the telescope position angles 
for the two 160 /im images differ by just 6° and so the PSF 
for the combined image is still tri-lobed). However, prelimi- 
nary modelling showed that an asymmetry of this magnitude 
could not be caused by the beam shape, and we concluded 
that this is a real feature that originates either in a back- 
ground source, or in a cold spot in the disk. The presence of 
such an extra component explains why the elliptical Gaus- 
sian fitting resulted in a higher inferred inclination from the 
160 /j,m image (26°) than at shorter wavelengths (20°). 

2 A. 2 SPIRE 

Follow-up observations were obtained on 2011 January 9 
with SPIRE (Spectral & Photometric Imaging REceiver, 
Griffin et al. 2010) using small-map mode, resulting in simul- 
taneous 250, 350 and 500 /xm images (see Table 1). The data 
were reduced using HIPE (version 7.0 build 1931), adopting 
the natural pixel scale of 6, 10, 14 arcsec at 250, 350 and 
500 /im respectively; the beam size at these wavelengths is 
18.2, 24.9, 36.3 arcsec. The resulting images are shown in 
Fig. 1. 

Emission is confidently detected at 250 and 350 /xm in 
the vicinity of the star, and is also found to be extended 
beyond the beam at these wavelengths. Surprisingly though, 
the emission is extended along a different position angle to 
the extension seen by PACS, and furthermore the peak is 
found to be offset significantly (~ 7 arcsec) from the star 
(see Fig. 1). For example, a 2D Gaussian fit to the 250 and 
350 /xm images finds FWHM for the major and minor axes 
and position angles and centres of: 43", 22", 82°, at (5.7"E, 
4.9" N), and 48", 28", 85°, at (6.9" E, 5.3" N). 

Since the pointing accuracy of Herschel is still 2 arc- 
sec for the SPIRE observations, such a large offset is highly 
unlikely to be caused by a pointing error. Rather we prefer 
to attribute the observed morphology, both the offset cen- 
tre and the different position angle, to the superposition of 
emission from at least two sources. That is, the circumstel- 
lar component seen at all PACS wavelengths and the back- 
ground source 27" W of 61 Vir seen most prominently at 
160 /xm with a spectrum rising toward longer wavelengths, 
and probably also a source to the N of the star. This super- 
position is illustrated in Fig. (3, left). 

This superposition makes it problematic to derive fluxes 
for the circumstellar component at SPIRE wavelengths. In 
this section we present total fluxes measured within aper- 
tures large enough to encompass the emission from both 
sources, noting that a significant fraction of this is expected 
to arise from the background object (s). Within a 40 arc- 
sec radius aperture, the emission at 250, 350 and 500 /xm is 
120 =b 21 mJy, 55 =b 17 mJy, and 20 =b 12 mJy. Again, these 
uncertainties have to be enhanced by a 7% repeatability 
and calibration uncertainty (see SPIRE observer's manual). 



In §3 an alternative approach is used to separate the flux 
from the different components by modelling the image as a 
circumstellar disk with offset point sources. 

2.2 Ancillary data 

2.2.1 VLA 

To assess the possibility that there is extragalactic emission 
in the vicinity of 61 Vir we searched the literature for previ- 
ous radio observations. VLA observations on 1995 January 
15 taken as part of the NRAO VLA Sky Survey (Condon 
et al. 1998) detected a 50.9 =b 1.9 mJy source at 1.4 GHz 
at a location (13h 18m 23.22s, -18d 18m 51.9s), with posi- 
tional uncertainty of 1 arcsec. This source was noted to 
be extended beyond the 45arcsec beam, with a fitted (de- 
convolved) FWHM of 33.2 zb 2.0 arcsec along a major axis 
at a position angle of 3.7 =b 2.4° (compared with a FWHM 
of < 18.4 arcsec along the minor axis). 

A search of the VLA archive found a 5 GHz scan of 
the region dating back to 1987 February 24 with ~ 10 arc- 
sec resolution (in DnC configuration). These data, imaged 
using uniform weighting (robust=-5), are shown in Fig. (3, 
middle) indicating a double-lobed structure. Fitting this as 
two Gaussians gives 8.0 mJy for the southerly lobe (13h 
18m 23.14s, -18d 19m 04.7s) and 8.8 mJy for the northerly 
lobe (13h 18m 23.13s, -18d 18m 39.3s), both of which have 
deconvolved sizes of around 7x4 arcsec; similar results are 
obtained from images reduced with natural weighting. The 
peak of the 1.4 GHz emission lies, to within the pointing 
uncertainty, directly between the lobes. 

At the epoch of the first (5 GHz) radio observations 
the star was offset ^ 39 arcsec NE of the centre of the radio 
emission. The high proper motion of 61 Vir, —1.07 arcsec/yr 
in RA and —1.06 arcsec/yr in Dec, means that the star had 
moved 11.9 arcsec to the SW by the time of the second 
(1.4 GHz) radio observations. Thus the lack of change in 
the centroid of the radio emission from 1987 to 1995 implies 
that the radio emission was not co-moving with the star 
(see Fig. 3); it can thus be assumed to be a background 
object. Indeed, the morphology and spectral index of the 
radio emission are typical for a non-thermal jet, which most 
likely originates from an AGN (Condon et al. 1998) that is 
directly between the lobes. 

The stellar proper motion has the unfortunate conse- 
quence that the star's location at the epoch of the Herschel 
observations (see §2.1.1) was just ^ 4 arcsec E of the puta- 
tive AGN. Thus it must be considered whether the emission 
seen with Herschel is associated with the extragalactic radio 
emission rather than being circumstellar. It is not possible 
to predict the far-IR brightness of the extragalactic emission 
from its radio flux, since the far-IR/radio correlation does 
not hold for radio- loud AGN. 

2.2.2 MIPS 

The possibility that the emission interpreted as circumstel- 
lar is in fact a background object is ruled out by archival 
MIPS observations taken at both 24 and 70 /xm on 2004 June 
21, that we have re-reduced here using an updated pipeline. 
Since the star contributes the majority of the 24 /xm emis- 
sion, this image does not provide information on the mor- 
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Figure 3. Disentangling the confusion toward 61 Vir. (Left) Superposition of the 250 jim contours onto the 160 fim image suggests that 
the 250 jim morphology is comprised of a circumstellar component (centred on the asterisk) and two offset sources (centred on the white 
dots). (Middle) The high proper motion of 61 Vir (see asterisks for the location of the star at the different observed epochs) confirms 
that double-lobed radio emission that is centred just 5 arcsec from 61 Vir at the epoch of the Herschel observations is background and 
not comoving with the star. The presented image shows 5 GHz observations from 1987, and the brown diamond shows the centroid of an 
extended 1.4 GHz source detected in 1995. (Right) The high proper motion of 61 Vir confirms that the extension of the 70 fim emission 
seen in 2011 with PACS (image reproduced from Fig. 1) was also seen in the archival Spitzer MIPS 70 jim image from 2004, which is 
shown with yellow contours at 50%, 75% and 90% of the peak value at 6.84 mjy/pixel. The high proper motion of the star confirms 
that this extension is comoving with the star, and so validates the circumstellar interpretation for the majority of the emission seen at 
70-160 /im. 



phology of any circumstellar emission. However, the 70 /xm 
image can be compared directly with that obtained at the 
same wavelength with Herschel 7 years later. 

The MIPS 70 /im observation has lower resolution 
(FWHM of 19 arcsec) than that of PACS (FWHM of 5.6 arc- 
sec), but the MIPS contours (e.g., those shown in Fig. 3 
right) are significantly extended beyond the beam, a point 
already noted in Lawler et al. (2009). Both 70 /xm images 
include the flux from the point-like star, though at 48.5 mJy 
this contributes only around 20% of the total emission, and 
the morphology seen by MIPS does not readily identify this 
as a distinct component as does the higher resolution PACS 
image. However, the MIPS 70 /xm emission is seen to be ex- 
tended along a position angle that is very similar to that 
seen in the PACS images, i.e., 56 =b 5°, with a 2D Gaussian 
resulting in FWHM of 24.5 =b 0.6" and 20.8 ± 0.5" along the 
major and minor axes, respectively. 

As previously noted the Herschel pointing accuracy is 
2 arcsec; to allow a proper comparison, the astrometry of the 
PACS image in Fig. (3, right) is that given by the Herschel 
pointing (i.e., the small correction used to recentre the image 
in Fig. 1 has not been applied). The pointing accuracy of 
Spitzer ^ is 1 arcsec, and again this has not been corrected 
for in Fig. (3, right), explaining the 2.7 arcsec offset of the 
peak from the stellar location at this epoch. We conclude 
that the stellar proper motion allows us to confirm that the 
excess emission is centred on the star, and that both epochs 
see the same emission morphology that is comoving with 
the star; i.e., it is not possible to interpret this emission as 
a stationary background object. 



^ http://irsa.ipac.caltech.edu/data/SPITZER/docs/ 
spitzermission /missionoverview / spitzertelescopehandbook / 12/ 



Given the lower resolution of the Spitzer beam, the 
MIPS 70 /xm image is not used in subsequent modelling of 
the dust distribution, although the total flux derived from 
the image using a 36 arcsec radius aperture was reported in 
§2.1.1 as 237 ±10 mJy, where it was compared with that 
derived from PACS. Virtually all the flux at this wavelength 
is within 20 arcsec radius, since smaller aperture radii mea- 
sure 223m Jy (16 arcsec radius) and 233m Jy in (20 arcsec 
radius) . 

Photometry on the 24 /xm image was used in the SED 
modelling. The total flux derived from fitting a PSF to the 
24 /xm image is 457.8 ± 4.6 mJy, compared with a stellar 
flux of 441.4 ±8 mJy, with no evidence for extension beyond 
the 5.5 arcsec beam. Thus this re-reduction has resulted in 
a modest upwards revision of the flux from the 450 mJy 
reported previously (Bryden et al. 2006; Trilling et al. 2008), 
and consequently there is now marginal evidence, taking un- 
certainties in the photospheric contribution into account, for 
a small ~ 1.8cr excess at this wavelength. This flux is consis- 
tent with that seen at 24 /xm by IRS (see §3.2 and Fig. 5). 
There is no evidence for any flux at the location of the pu- 
tative AGN, which was 13 arcsec from the star at the time 
of the observations; isolated point sources above 0.51 mJy 
could have been detected in the MIPS image with > 3cr con- 
fidence, but given its proximity to the star the limit on the 
24 /xm AGN flux is ^ 5 times higher. 

2.2.3 STIS 

Archival HST/STIS images were searched for evidence ei- 
ther for scattered light from the debris disk or for extra- 
galactic emission. The data were obtained as part of a cali- 
bration program (CAL/STIS-8419, PI Charles Proffitt), and 
were not optimized for a debris disk search. For example. 



images were obtained with a 1100 x 588 pixel subarray 
(5.58^^ X 2.98^^) and contemporaneous images for the subtrac- 
tion of the point spread function were not obtained. 61 Vir 
was observed on 2000 March 11 and again on 2000 May 14. 
In both observations 61 Vir was occulted at two positions 
behind a tapered focal plane occulting wedge. At the first 
position the wedge is 1.0^^ wide, and 160 frames with 7 sec- 
onds integration each were obtained. At the second position 
the wedge is 1.8^^ wide, and 56 frames of 7 seconds each were 
acquired. The only difference between the two epochs of ob- 
servation, aside from wedge width and total exposure time, 
is that in 2000 March 11 the telescope was oriented such that 
north is 102.44° counterclockwise from the +Y axis of the 
detector, whereas in 2000 May 14 north is 51.19° clockwise 
from the +Y axis. For PSF subtraction we registered the two 
epochs of data (for a given wedge position) and subtracted 
61 Vir from itself. The difference in telescope position an- 
gle means that a linear feature symmetric about the star, 
such as a debris disk midplane, would be offset from itself 
by 26.4°. 

Fig. 4 shows our PSF-subtracted image with 61 Vir oc- 
culted behind the 1.8^^ wide occulting wedge position. Neb- 
ulosity from circumstellar dust is not detected. The obser- 
vations made with the 1.0^^ wide occulting wedge position 
are more sensitive closer to the star, and from these we esti- 
mate an approximate 3a sensitivity for an extended source 
at 3.7" radius of 20.0 mag arcsec"^ in the VEGAMAG V- 
band. At radii > 10^^ the data for both occulting wedges have 
an approximate 3a limiting sensitivity of 21.5 mag arcsec"^. 
When combined with the far-IR images, these limits set use- 
ful constraints on the albedos of the dust grains in the disk 
that are discussed in §3.3.1, 3.3.2 and 3.3.3. Unfortunately, 
given the epoch of the STIS observations, the field of view 
does not encompass the possible extragalactic source that is 
currently west of 61 Vir. 



3 MODELLING 

Here we consider the structure of the disk required to explain 
both the images presented in §2, and the spectrum of the 
disk emission. To do so, first the spectrum of the star is 
determined (§3.1), then additional data needed to constrain 
the shape of the disk's emission spectrum is described (§3.2), 
before modelling of all the observational data is undertaken 
in §3.3. 



3.1 The star 

The main sequence G5 star 61 Vir has an age estimated as 
4.6 ± 0.9 Gyr using gyrochronology (Wright et al. 2011; Vi- 
can 2012), a mass of O.88M0 and a metallicity of [Fe/H] = 
-0.02 zb 0.01 (Sousa et al. 2008). Its spectrum was de- 
termined by fitting PHOENIX models from the Gaia grid 
(Brott & Hauschildt 2005; Hauschildt k Baron 2009) to 
the observational data available for the star at wavelengths 
shorter than 12 /xm where we can be confident that there is 
minimal contribution from the disk. This fit includes fluxes 
from Tycho2 (Bt,Vt), 2MASS (J,H,Ks), AKARI (9 /xm), 
IRAS (12 /im) (see Kennedy et al. 2012 for more details on 
this procedure). 
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Figure 4. HST/STIS optical image of 61 Vir. North is up, east is 
left, and the white circle indicates 10^^ radius. We show the data 
with a 3 X 3 pixel median filter to improve sensitivity to faint 
extended sources. A software mask blocks orthogonal diffraction 
spikes and the occulting bar that differs from a horizontal orien- 
tation by 12°. With a position angle of 65°, the midplane of the 
61 Vir debris disk is offset from the bar by another 13 degrees 
clockwise, but is not detected in this image. 



A best fit was found with Tgff = 5602K scaled to a lumi- 
nosity of O.84L0. Photospheric fluxes based on this model 
have a la uncertainty of ±2%. The same process was ap- 
plied to other stars from the DEBRIS survey (Kennedy et 
al. 2012), and a comparison of the predicted fluxes to those 
observed at longer wavelengths for stars where a disk is not 
present provides confidence that this uncertainty is reason- 
able. The predicted stellar spectrum is shown in blue on 
Fig. 5, and fluxes predicted in the Herschel bands, with no 
colour correction, are 50.9, 25.1, 10.0, 4.0, 2.0, 1.0 mJy at 
70, 100, 160, 250, 350 and 500 /xm, respectively, while for 
MIPS 24 /im the stellar flux is predicted at 441.4 mJy. 

3.2 Excess spectrum 

The spectrum of the excess (i.e., the observed flux less that 
expected from the photosphere) is plotted in Fig. 5. In ad- 
dition to the fluxes from §2, this plot includes the IRAS 
60 /im flux measured using SCANPI^ of 210 ± 60 mJy, and 
importantly also includes the Spitzer IRS spectrum (AOR 
15998976) that was first presented in Lawler et al. (2009). 
Here we use the more recent CASSIS reduction (The Cor- 
nell AtlaS of Spitzer/IRS Sources, Lebouteiller et al. 2011), 
except that the spectra from the individual modules were 
aligned and then the whole spectrum was scaled so that 
the flux in the 8-9 fim range matched that of our model 
stellar spectrum (§3.1). This is because the 8% calibration 
accuracy of IRS (Lawler et al. 2009) is lower than that of 
our photospheric model, and the shape of the IRS spec- 
trum indicates that there is no contribution from the disk 
at such wavelengths (e.g., Lawler et al. 2009 report an ex- 
cess of —3.6 ± 1.4mJy from the 8.5-12 fim range). To aid 

^ http:/ /scanpiops.ipac.caltech.edu:9000 
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wavelength (|im) 

Figure 5. Spectral energy distribution of 61 Vir. Total fluxes 
are black circles, open circles are predicted photospheric fluxes, 
grey circles are photosphere subtracted (i.e., excess) fluxes; upper 
limits are plotted as downward pointing triangles (grey and open 
for excess fluxes and black and solid for total fluxes). Data points 
from the pre-binned IRS spectrum are shown in green beyond 
10 fim. The photospheric model is shown with a blue line, and is 
a flt to data < 12 fim.. The longer wavelength excess fluxes are 
fltted as a 68K black body of fractional luminosity 2.7 x 10~^ that 
is shown in red. The total flux from the models (photosphere plus 
black body) is shown with a black solid line. 



the modelling, the IRS spectrum was averaged over 7 bins 
2 — 4 /xm wide, with the uncertainty in the measurement 
taken from the distribution of measurements within that bin 
summed quadratically with the 2% calibration uncertainty 
of the spectrum (Beichman et al. 2006). For example, this 
process finds a flux of 270 =b 10 mJy in the range 30-34 /xm, 
which is 24.4 mJy higher than the predicted photospheric 
flux, and can be compared with a 30.8 ± 2.0 mJy excess 
reported in Lawler et al. (2009). 

The overall shape of the excess spectrum resembles that 
of a black body. Indeed the whole spectrum can be fitted 
reasonably well by a single temperature at 68K with a frac- 
tional luminosity / = 2.7 x 10~^ (see Fig. 5). This fit also 
highlights where the spectrum departs from a black body. 
The most prominent discrepancy occurs at > 160 /xm, where 
the observed fluxes are higher than those of the black body 
model. While such discrepancies could indicate that the disk 
spectrum contains multiple temperature components, or in- 
deed features associated with specific dust compositions, it 
is more likely that these longest wavelength fluxes are con- 
taminated with what is assumed to be background emission, 
since they include all flux within 40 arcsec of the star (see 
§2.1.2); these contaminating sources are accounted for in the 
modelling of §3.3. 

One thing that is clear from the SED, however, is that 
there is no evidence for a hot (^ 68K) component suggesting 
that the inner regions are relatively empty. This is supported 
by the star subtracted 70/im image (Fig. 2). Although an 
inner hole is not evident as a decrease in surface brightness 
in the central pixels, the resolution is such that the central 
pixels are expected to include flux from the outer regions of 
the disk. However, the image would be expected to be much 



more centrally peaked if there was dust close to the star, 
and this point is quantified through modelling in §3.3. 

3.3 Models 

The aim of the modelling is to arrive at an axisymmetric 
circumstellar disk structure that fits all of the images in 
Fig. 1 as well as the spectrum in Fig. 5. Clearly this will not 
be possible at the longer wavelengths without also includ- 
ing background sources. Given the interpretation in previ- 
ous sections, the simplest morphology for the background 
sources is that there are two point sources, one ^ 27 arcsec 
to the west, and another ^10 arcsec to the north. Thus, in 
the models all images include such sources at offsets from 
the star that are free parameters (but that are the same 
for all wavelengths) , with fluxes at each wavelength that are 
also free. All of these parameters are determined by a best 
fit to the images, however we concluded that the quality of 
the 500 /im observations was not sufficient to derive fluxes 
for the individual background sources at this wavelength. 
The models also include emission from the point-like star at 
the level derived for the photosphere (§3.1). 

Different models are tried for the disk profile, discussed 
in more detail below. However, there are some common ele- 
ments between all models, and how the models are assessed, 
that are described here. First, no attempt is made to con- 
strain the grain properties or size distribution, since in the 
absence of spectral features this will be degenerate with 
other assumptions. Rather it is assumed that the emission 
spectrum of the disk at any given radius resembles a mod- 
ified black body that is defined by its temperature, T, and 
the wavelength beyond which the spectrum is modified, Ao, 
and the modification to the spectral slope f5 (i.e., the spec- 
trum has an additional factor oc (Ao/A)^ for wavelengths 
A > Ao). The parameters Ao and /3 are left free, but given 
the level of confusion at sub-mm wavelengths they are poorly 
constrained, as discussed further below. The temperature is 
assumed to be T = /xThb, where Tbb = 266r~^^^ is the tem- 
perature (in K) that black bodies would attain at a distance 
r (in AU) for this stellar luminosity, and /t is a free param- 
eter. It was not found to be necessary to change the —1/2 
power law radial dependence of temperature in the models. 
In models that include multiple components, the different 
components were sometimes allowed to have different prop- 
erties (/t, Ao, P). 

The emission from different annuli in the disk model was 
then determined by the vertical geometrical optical depth 
distribution T(r), that had different assumptions in the dif- 
ferent models. An image of the model was made by integrat- 
ing the flux along the line-of-sight of different pixels given 
2 free parameters that define the orientation of the orbital 
plane (the inclination to our line-of-sight, /, and the position 
angle, PA); if the pixels undersampled the disk structure a 
denser grid of pixels was used, and then rebinned, to en- 
sure that all of the disk flux was captured. The background 
sources and the star were added to individual pixels, and 
the whole image convolved with an appropriately oriented 
telescope PSF (e.g., Kennedy et al. 2012). 

To quantify the success of the model, the model images 
were subtracted from the observed images, then divided by 
the la pixel-to-pixel uncertainty given in Fig. 1, following 
which the resulting image was squared. This was then used 
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to derive a by summing over pixels near the disk^. A 
was also derived from the SED, but only including the 
wavelengths where the photometry was not derived from one 
of the Herschel images, i.e., including the AKARI 18 /xm, 
MIPS 24 /xm, and the IRS spectrum (but not the MIPS 
70 /xm as previously noted). 

The for different images and that from the SED 
were then summed using the following weighting. Each pho- 
tometric data point from the literature received a relative 
weight of 1. Images were weighted by the area above the 3a 
flux level divided by the area above the same flux level of 
the PSF scaled to the image peak. With this prescription, 
an unresolved image has the same weight as a photometric 
data point, and resolved images are weighted by the number 
of resolution elements that the source covers. 

Since the location of the star in the images is not well 
constrained, this location is also left as a free parameter for 
the 3 observations (100/160, 70/160, and 250/350/500). The 
free parameters of the model are thus: 6 for the location of 
the star, 4 for the offsets of the two background sources, 
12 for the fluxes of the background sources at the 6 wave- 
lengths, 2 for the orientation of the disk plane, 3 for the 
dust properties (/t, Aq, /3), plus those describing the spatial 
distribution. 

With so many parameters, finding the best fit by grid 
methods is prohibitive. Instead, we initially find reasonably 
fitting models by varying parameters by hand, which are 
further refined with least squares minimisation. The min- 
imisation step means that our models at least lie at a local 
minimum in parameter space, but does not exclude the pos- 
sibility of a better fit in a different region of parameter space. 
Nevertheless, since the models provide a reasonable fit to the 
observations, we believe that they are at least representative 
of the variety of models that can fit the data. However, we 
prefer not to trust the uncertainties in the model parame- 
ters returned by the minimisation process, and quote these 
based on other considerations as outlined below. 



3.3.1 Single extended component 

The simplest model tried had a radial profile given by r (r) = 
T"o(r/rmin)^ bctwccn inner and outer radii rmin and rmax- 
The best fit found was for rmin = 29 AU, rmax = 350 AU, 
p = -1.1, TO = 6.9 X 10-^ / = 77°, /t = 1.9, Ao = 220 /xm, 
f5 — 1. The background sources were found to be offset by 
(-27", +2") and (-4", +7") with fluxes resembling cool emis- 
sion with temperatures of < 40K (see Fig. 6), and so are 
consistent with extragalactic emission. The resulting model 
images are shown in Fig. 7, that can be compared directly 
with Fig. 1, since the figures employ the same colour scale 
and contours. Fig. 7 also shows the residual image (i.e., ob- 
served minus model). Clearly the model provides an excel- 
lent fit to the observations, and the residuals in the vicinity 
of the star are consistent with the noise elsewhere on the 
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Figure 6. Model SEDs of 61 Vir. The top panel shows the con- 
tinuous 29-350AU dust distribution model (§3.3.1), the middle 
panel shows the 40 and 90AU ring model (§3.3.2), and the bottom 
panel shows the collisionally evolved pre-stirred 1-350AU model 
(§3.3.3). Fluxes shown in black are those observed less that of 
the photosphere, i.e., these are the same fluxes plotted in grey on 
Fig. 5, except that here detections that are below 3cr are plotted 
both at their measured value (with their Icr uncertainty) and as 
triangles to indicate the 3cr upper limit. The model disk fluxes are 
shown with a thick red line, and contributions of the two com- 
ponents of the bottom two models have also been separated out 
with dotted and dashed lines, for the innermost and outermost 
component respectively. The inferred spectra of the background 
sources are shown in blue. 



^ Near here means those pixels likely to contain non-zero flux. 
These were chosen by taking an observed image and setting all 
pixels over 3cr to 1, and setting those below 3cr to 0, then con- 
volving with the PSF and retaining pixels > 0.1. This procedure 
yields a mask that is somewhat larger than the disk, ensuring 
that models are also penalised for being too extended. 



image; there is no evidence for any asymmetry in the disk 
structure. The resulting reduced was 1.51. 

Although the algorithm provides uncertainties on the 
model parameters, the most instructive way to qualify how 
well the radial distribution is constrained is to compare this 
result with an alternative model that provides a compara- 
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Figure 7. Model images of the PACS 70 iim (left), 100 iim (middle) and 160 /im (right) emission distribution toward 61 Vir using 
the single extended component model of §3.3.1 which has a continuous dust distribution from 29-350AU. The top row of images can be 
compared directly with Fig. 1, since these are model observations, and employ the same contours and colour scale. The bottom row are 
the residual images (observed minus model), with a colour scale that goes from —2 to +6 times the Icr uncertainty in each observation, 
and contours are shown in black at +4cr. The white contour shows the region where the observed flux is above 4cr (from Fig. 1). 



ble fit to the image (§3.3.2). However, both models come 
to similar conclusions on the parameters defining the dust 
properties. These show that the factor by which the dust 
temperature exceeds that of black bodies in thermal equilib- 
rium, /t, is well constrained, since it is essentially set by the 
emission spectrum. This factor is in line with expectations 
that it should be greater than unity, because the ^ 1 /xm 
dust grains that dominate the cross-sectional area of a de- 
bris disk absorb starlight efficiently but reemit inefficiently 
at longer wavelengths (e.g.. Fig. 3 of Bonsor & Wyatt 2010). 
It is also in line with observations of other disks that have 
been imaged (e.g.. Fig. 9 of Rodriguez & Zuckerman 2012), 
since one would expect disks to be around times larger 
than would have been anticipated from their temperature 
assuming the grains to act like black bodies. On the other 
hand the factors Ao and f5 are poorly constrained, since they 
depend on the fraction of the emission seen in the 350 and 
500 /im images that comes from the disk, and the resolu- 
tion is insufficient to separate the disk emission from the 
background sources (particularly that to the north), which 
moreover may have a more complicated structure than the 
2 point source model that has been assumed; the circum- 
stellar emission will be less confused in > 5 — 10 years once 



the star has moved away from the background sources (see 
Fig. 3 middle). The uncertainty on the disk inclination is 
estimated as ±4°. 



This model was also used to predict the correspond- 
ing scattered light images, under the assumption that the 
grains scatter light isotropically and have a constant albedo 
uj at all optical wavelengths at which they interact with 
starlight. These model images were then compared with 
the constraints on the optical surface brightness set by the 
STIS images in §2.2.3 to set constraints on uj. The scattered 
light in this model is brightest at its inner edge at around 
3.7arcsec, constraining the albedo to uj < 0.31. This is not 
a definitive statement about the grain albedo, since other 
disk structures also fit the thermal emission and result in 
slightly different constraints on the albedos (see §3.3.2 and 
3.3.3). However, we note that u < 0.31 is consistent with 
the grain albedos inferred for other debris disks imaged in 
scattered light around solar type stars (e.g., Backman et al. 
2009, Krist et al. 2010, Golimowski et al. 2011), as well as 
for Solar System interplanetary dust particles (Kelsall et al. 
1998). 
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Figure 8. Residual images (observed minus model) of the PACS 70 /im (left), 100 /im (middle) and 160 /im (right) emission distribution 
toward 61 Vir. The top row shows the residuals of the best fit model of §3.3.2 which has two narrow rings at 40 and 90AU. The bottom 
row shows the residuals of best fit model of §3.3.3 which corresponds to a collisionally evolved pre-stirred disk that initially extended 
from 1-350AU. The colour scales and contours are the same as Fig. 7. 



3.3.2 Two radii 

The next model had a radial profile given by two narrow 
rings at radii of ri, r2, each of which had a constant optical 
depth n and T2 across a width of lOAU, and with different 
particle properties /ti, Xoi, ft. The best fit found had ri = 
40AU, r2 = 90AU, n = 1.4 x 10"^, = 1.7 x 10"^, with 
particle properties /ti = 1-9, Aoi = 70 /xm, fii = 0.53, 
and /t2 = 1-5, A02 = 280 /xm, ft = 0.95, with background 
sources with fluxes and locations very similar to those of 
§3.3.1 (see Fig. 6). Again the model provides an excellent 
fit to the observations (see Fig. 8), with residuals that look 
broadly similar to those of Fig. 7, and has a total reduced 
of 1.94. 

As this model provides a comparable fit to the same set 
of observations, it illustrates the uncertainty in the radial 
profile. For example, we should not take too seriously the 
conclusion of §3.3.1 that the radial profile extends out as far 
as 350 AU, since the contribution from the outermost regions 
beyond lOOAU is negligible, even if the fit of a single compo- 
nent model is improved by including emission in these outer 
regions. The larger outer extent is favoured because there is 
emission out to 20" to the NE at 70-160 /xm. Fig. 1 shows 
this is not equally matched to the SW. Since the two ring 



model does not include emission at such large distances, this 
results in significant (> 4a) residual emission being present 
- 20 arcsec ENE of the star in Fig. 8. 

The resolution of the images is also not such that we 
can tell if the distribution between 30 and lOOAU is con- 
tinuous or concentrated in narrower features. Similarly the 
inner edge is not well known, but to explore this we consider 
a third model in §3.3.3. We do not attribute any physical sig- 
nificance to the fact that the two rings are inferred to have 
different Ao and /3, which arises to ensure an appropriate 
balance of flux from the two components at all wavelengths. 
Rather, again, we consider these parameters to be poorly 
constrained. With this model, the non-detection in the STIS 
images result in a constraint on the dust albedo of cj < 0.25. 
This is sightly more stringent than the model of §3.3.1 as 
the emission in this model is more concentrated. 

3.3.3 Collisionally evolved pre-stirred disk 

This model is the same as that of §3.3.1, except that we also 
allow the optical depth distribution to extend continuously 
inward of rmin to an inner edge at ro with power- law index po 
(which is positive so that optical depth increases with radius 
in this region). Thus, the disk surface density increases from 
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an inner edge (at ro) to a maximum (at rmin), then turns 
over and decreases towards the outer edge (at rmax)- 

Given that we have added parameters to a model that 
already had sufficient parameters to explain the observa- 
tions, we do not intend to present a full exploration of this 
parameter space. Rather, the intention here is to test a hy- 
pothesis that the planetesimal disk initially extended all the 
way from an inner edge just outside the known planets at 
ro (at around lAU) out to some large distance rmax, and 
that the inner regions of this disk (inside rmin) have since 
been depleted by collisional erosion, whereas the outer re- 
gions (beyond rmin) have yet to reach collisional equilibrium. 
If the disk is pre-stirred (see Wyatt 2008), in the sense that 
collisions between the planetesimals were destructive from 
shortly after the star reached the main sequence, then the 
inner regions would be expected to attain an optical depth 
profile with po = 7/3 (Wyatt 2008; Kennedy & Wyatt 2010; 
see §4.4). It is this model that we consider here. The physi- 
cal interpretation of derived parameters such as rmin is dis- 
cussed in §4.4, and a discussion of the self-stirred (Kenyon & 
Bromley 2010; Kennedy & Wyatt 2010) and planet- stirred 
(Mustill & Wyatt 2009) variants of this model is given in 
§4.5. 

This model has a reduced of 1.60, and as with the 
previous models we find an excellent fit to the data; e.g.. 
Fig. 8 shows residuals that are very similar to those of the 
single extended component model of Fig. 7. The inner edge 
ro is not well constrained; it can extend all the way to 
the star because the surface density of the (collisionally de- 
peleted) inner region increases very steeply with radius so 
little emission originates from the innermost part of the disk. 
The model we present here has ro = 1 AU at which radius 
the optical depth is tq = 7.4 x 10~^. The turnover occurs at 
^min = 43 AU, so has moved outwards compared to the in- 
ner edge of the single component model of §3.3.1 to account 
for the extra inner emission. The outer edge lies again at 
Tmax = 350 AU (but as noted above could be much closer). 
Other parameters are similar to the single component model, 
e.g., Ao = 225 /im, /3 = 0.7, /t = 1.8, largely because the 
model assumed an optical depth profile in the inner regions 
that was steep enough to mean that this component is not 
a major contributor to the disk emission (Fig. 6). 

Thus, the observations are consistent with a small level 
of emission being present within 30 AU. To quantify this fur- 
ther, we repeated this modelling, but decreased po, keeping 
ro at lAU, until the best fit with the assumed parameters 
was unacceptable. This procedure found that distributions 
with Po < 1 are ruled out, both because they result in a 
level of mid-IR emission that exceeds that seen in the IRS 
observations, and because the 70 /xm image becomes too 
centrally peaked. In other words, the slope of the inner edge 
is not well constrained, except that we can say that it rises 
sharply enough for there to be negligible fiux in the inner 
regions. From the IRS spectrum this means that the inner 
component must have a 10 fim flux below around 50 mJy 
(i.e., 2% of the stellar photosphere) and a 32 /xm flux below 
54 mJy assuming the outer component does not contribute 
any flux at this wavelength (which is unlikely so the true 
flux limit will be lower). From the 70 fim image, the inner 
component must contribute below around 6 mJy (the 3a 
point source uncertainty). The albedo constraint from the 
STIS non-detection is slightly higher than that of the previ- 



ous models at uj < 0.57, since the emission in this model is 
more diffuse. 



4 IMPLICATIONS FOR THE 61 VIR 
PLANETARY SYSTEM 

4.1 The 61 Vir planetary system 

The dataset of radial velocity measurements for 61 Vir has 
expanded since the observations reported in Vogt et al. 
(2010), allowing further consideration of the properties of its 
planetary system. Segransan et al. (in prep.) report 142 pre- 
cise HARPS radial velocities (averaged values per night) of 
61 Vir between JD = 53037 and 55948 that conflrm the ex- 
istence of the inner two planets, with parameters consistent 
with the values of Vogt et al. (2010) (i.e., a 5.IM0 planet 
at 0.05AU and a 18.2M© planet at 0.218AU). However, the 
outermost planet (i.e., the 22.9M0 planet at 0.478AU) has 
not been conflrmed in the HARPS data at this point. The 
temporal variation of the log {R'hk) activity index (which 
is also available from the HARPS spectra) shows that, al- 
though the star has been quiet most of the time, there was a 
burst of activity between JD =54800 and 55220. Such active 
periods can influence the interpretation of longer periodici- 
ties in the data, and so affect our ability to extract longer 
period planets. As the existence of this planet is not crucial 
to the conclusions of this paper, this issue is not discussed 
further here. 

Valuable constraints can also be set on the properties 
of planets that can be ruled out by the HARPS data (e.g.. 
Mayor et al. 2011). Since the 2-planet system fits the radial 
velocities to within K ^ 1 (in ms""*^), this crudely allows us 
to rule out planets on circular orbits at semimajor axes a^i 
(in AU) that are more massive than Mpi sini > lOK^a^i in 
M0. However the limit would be higher than this for plan- 
ets on periods much longer than the length over which the 
observations have been taken [r^ 8 years). Fig. 9 shows the 
limit derived from an analysis of the HARPS data (exclud- 
ing that obtained during the activity burst) from which the 
solution for the 2 inner planets has been subtracted. For each 
point in the mass-period diagram, a corresponding circular 
orbit is added to the residuals of the 2-planet model, with 
20 phase values (spaced equidistantly between and 2ti). 
For each phase, this process is repeated 1000 times with 
the residuals randomly permuted, and we consider that the 
planet is detected only if we see a signal in the resulting pe- 
riodogram that is higher than the 1% false alarm probability 
limit at all phases. The detection threshold obtained agrees 
well with that expected from the crude analysis above, and 
shows that the data exclude planets that are more massive 
than Saturn with semimajor axes inside 6AU. 

4.2 Interaction between known planets and debris 

Since the planet masses given in §4.1 were derived from 
radial velocities the masses reported are actually Mpisinz 
rather than Mpi, and the inclination of the planets' orbits 
to the sky plane (i) is unknown. Our observations determine 
the inclination of the disk's plane of symmetry to the sky 
plane to be 77°. If we assume that the planets' orbits are 
coplanar with the disk then their masses should be revised 
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Figure 9. Detection limit derived from the HARPS radial ve- 
locities. The blue line indicates the limit above which planets on 
circular orbits would have been detectable in the data (see text for 
details). The red dots indicate the position of the 2 inner planets. 



modestly upward by 3% from the values given above. Such 
a small change would not affect previous conclusions in the 
literature about the dynamics of this system (Vogt et al. 
2010; Greenberg & van Laerhoven 2012). 

However, the planets do not have to be coplanar with 
the disk, since the two are sufficiently spatially separated 
that there is little dynamical interaction between them. For 
example, although the planets' secular perturbations will 
make the orbits of disk particles precess so that their distri- 
bution would end up symmetrical about the invariable plane 
of the planetary system, the timescale for this to happen is 
at least the secular precession timescale. At 30 AU this is 
60 Gyr considering perturbations from the three Vogt et al. 
(2010) planets, and even longer if the outermost planet is 
excluded. This also precludes the disk having been stirred 
by these planets (Mustill & Wyatt 2009); i.e., the collisional 
cascade in the disk must have been ignited by some other 
mechanism (see §4.5). Rather, if there is to be an interac- 
tion between the known planets and the disk there would 
need to be planets in the intervening 1-30AU region, either 
to allow disk material to be scattered to these inner regions 
(e.g., Bonsor & Wyatt 2012), with implications for the hab- 
itability of the inner planets (e.g., Horner & Jones 2010), or 
for the planets' gravitational perturbations to influence the 
intervening planets which in turn influence the outer regions 
(e.g., Zakamska & Tremaine 2004). 

The assumption that the planets and the debris disk are 
coplanar is justified by the premise that the debris disk and 
planets were born from the same protoplanetary disk, which 
had a single plane of symmetry. Nevertheless this does not 
guarantee that the two should remain coplanar, since the de- 
bris disk may have been perturbed by passing stars, and the 
planets could have undergone a dynamical interaction that 
misaligned their orbital planes (Naoz et al. 2011). The spin 
axis of the star can also become misaligned with the proto- 
planetary disk under certain conditions (Lai et al. 2011). In 
this respect it is worth noting that the stellar rotation period 



of 29 days and vsini of 1.6 =b 0.5 km/s suggest the stellar 
pole is inclined 72° to our line-of-sight, i.e., consistent with 
being nearly perpendicular to the debris disk's symmetry 
plane. Thus there is no evidence at the moment to suggest 
that the angular momentum vectors of the star, planets and 
debris disk are misaligned. 

4.3 Intervening planets? 

Given that an axisymmetric disk model provides a good fit 
to the observations (§3), disk asymmetries cannot be in- 
voked as evidence for more distant planets as has been the 
case for other imaged disks (e.g., Wyatt 2008; Krivov 2010). 
However, the emptiness of the inner region does require an 
absence of planet esimals, and this emptiness can be read- 
ily explained if there are additional planets in the 1-30AU 
region, since such planets could dynamically deplete the re- 
gion (e.g., Lecar et al. 2001). In §4.1 we already used the 
radial velocity data to set constraints on the planets that 
are allowed to exist within 6AU (and have remained unde- 
tected, see Fig. 9). It is not possible to set hard constraints 
on more distant planets, except for generalities such as that 
the planets cannot reside too close to the inner edge of the 
disk (Quillen 2006; Mustill k Wyatt 2012). However, in §4.4 
we consider the possibility that no planets exist in the 1- 
30 AU region. 

4.4 Plausibility of pre-stirred extended disk model 

In §3.3.3 it was suggested that the observed depletion could 
arise from collisional erosion, i.e., without requiring any 
planets to exist in the 1-30AU region. To confirm that this 
is a plausible explanation, we show here that the observa- 
tionally inferred parameters for this model (rmin, tq and 
p) correspond to physically plausible parameters within the 
context of this scenario. 

First we consider the regions beyond rmin which, in this 
interpretation, have yet to deplete, so can be used to deter- 
mine the initial mass distribution of the debris disk; we will 
assume this initially extended all the way in to ro according 
to S = So(r/ro)^, where S is the mass surface density (e.g., 
m M® AU"2). The mass surface density simply scales with 
the optical depth according to the size distribution, which 
we assume to be a single power law with the —3.5 exponent 
expected (under certain ideal assumptions) for collisional 
equilibrium, from a maximum planetesimal size Dmax (in 
km) , down to the size at which dust is removed by radiation 
pressure. Using equation 15 of Wyatt (2008), for a O.88M0, 
O.84L0 star, assuming a solid density of 2700 kg m~^, 

S/r = 0.19DK^^ (1) 

in M® AU"^ 

For reference in later paragraphs, once we have esti- 
mated Dmax, we will compare the mass surface density de- 
rived with that of a fiducial minimum mass solar nebula 
(MMSN), which contains I.IM® AU~^ of solid material at 
lAU, with surface densities falling off to larger distance as 
(r/ro)"^/^ (e.g., Weidenschilling 1977). The best fit model 
of §3.3.3 has an optical depth at 43AU of 4.8 x 10~^ which is 
around an order of magnitude higher than that in the Kuiper 
belt (see Fig. 13 of Vitense et al. 2010), and corresponds to 
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a mass surface density at that location of 2.3 x 10 ^-Dmax 
times that in the MMSN. 

The plausibihty of this interpretation comes down to 
whether it is reasonable to expect the transition to colli- 
sional equilibrium to occur at 43AU given the stellar age of 
4.6 ±0.9 Gyr. As well as the aforementioned parameters, the 
collisional lifetime of the largest planetesimals in the disk de- 
pends on their strength and eccentricity e (which sets 
collision velocities). For the assumptions above, this colli- 
sional lifetime can be estimated (see eq. 13 of Wyatt et al. 
2007a) to be 

tc{D^,.) = 1.4 X 10-'r'^'[Qh'^'DK^^e-'^']T-' (2) 

in years. In the inner regions that are assumed to be in 
collisional equilibrium this lifetime should be the age of the 
star, explaining why we set po = 7/3 in §3.3.3. 

Thus, for the collisional lifetime at 43AU to be equal to 
the 4.6Gyr age of the star, the quantity in square brackets 
in eq. 2 must be of order 2.4 x 10^ (±50% due to the uncer- 
tainties in stellar age and disk optical depth and turnover 
radius). This quantity has been estimated from the detec- 
tion statistics as a function of age for the ensemble of debris 
disks around A stars to be 7.4 x 10^, with a dispersion of 
Idex about this average disk value (Wyatt et al. 2007b). 
For sun-like stars the same process gives this quantity as 
2.9 X 10^ with a similar dispersion (Kains, Wyatt & Greaves 
2011). So it appears that, for this interpretation to be valid, 
the properties of planetesimals around 61 Vir would have to 
be broadly similar to those of average A star debris disks, 
but would be in the 2a minority when compared with debris 
disks known around sun-like stars. The sense of this differ- 
ence is that the planetesimals would have to be smaller, or 
on more elliptical orbits. We find that the required quan- 
tity is obtained for Dmax = 2 — 5 km for e = 0.05 — 0.1 
for two different strength laws, i.e., that of weak aggregates 
(Stewart & Leinhardt 2009) and that for low temperature 
ice (e.g., Vitense et al. 2010). 

To conclude on the pre-stirred model, we find that the 
observations can be explained if the disk is made up of plan- 
etesimals up to a few km in size, with a mass surface density 
around 0.4% that of the MMSN. Despite its low mass, which 
is comparable with that of the Kuiper belt, the disk in this 
model is detectable because the mass is locked up in rela- 
tively small planetesimals. Such a low mass is somewhat at 
odds with the presence of planets, since if these have even 
modest rocky cores, then the solid mass surface density in 
the inner regions is at least an order of magnitude above the 
MMSN. Even if all of the mass inside 43AU in the hypo- 
thetical initial debris disk was placed in the planet region, 
that would only correspond to O.2M0. One would therefore 
have to conclude that the disk was initially much more mas- 
sive, and that the majority of the mass has been depleted 
somehow, such as by inwards migration due to gas drag (or 
simply never having grown into planetesimals). 

Remember though, that this conclusion is based on the 
premise that the 43AU peak arises from collisional evolution. 
If the inner edge is truncated by some other mechanism, then 
the debris disk could include much larger planetesimals and 
so be much more massive. For example, if the planetesimals 
were up to 2000km in diameter the debris disk's mass surface 
density would be inferred to be around 10% of a MMSN at 
43AU. Given the uncertainties, this initial mass is more com- 



patible with the known planets, but would still require the 
original protoplanetary disk to be significantly more mas- 
sive than the remnant debris disk, and likely also that the 
planets we know today include a reasonable fraction of the 
mass that was originally inside 43 AU. 

4.5 Disk stirring mechanism 

In the model discussed in §4.4, the disk was assumed to 
be stirred from early on, since that allowed us to equate 
the collisional lifetime (eq. 2) with the stellar age. Such a 
model is called pre-stirred in Wyatt (2008), and makes no 
assumptions about the stirring mechanism. 

A variant on the pre-stirred model, that includes a self- 
consistent consideration of the stirring mechanism, is the 
self-stirred model. In this model collisions in an extended 
planetesimal belt initially occur at low enough velocity for 
planet esimal growth to occur (Kenyon & Bromley 2010; 
Kennedy & Wyatt 2010). Such models allow planet forma- 
tion processes to continue throughout the several Gyr life- 
time of the star. Within such a model a peak in surface 
brightness naturally corresponds to the location at which 
planetesimals have recently grown to around Pluto-size, 
since at this time gravitational scattering stirs neighbour- 
ing planetesimals resulting in enhanced dust production. For 
the growth of Pluto-size objects, i.e. those large enough to 
stir the disk, to have reached 43 AU after 4.6Gyr of evolu- 
tion, this would require the original protoplanetary disk of 
61 Vir to be less than 0.01% of a MMSN (see equation 9 of 
Kennedy & Wyatt 2010). Since this would require the disk to 
be made up of entirely of sub-km-sized planetesimals, which 
is contrary to the premise of the model that planetesimals 
are growing up to Pluto-size, we do not consider this to be 
a viable explanation. 

Fully formed planets that are well separated from the 
disk can also stir the disk through their secular perturba- 
tions, i.e., the planet- stirred model (Mustill & Wyatt 2009). 
We already showed in §4.1 that the known planets cannot 
have stirred the disk given the age of the star. However, 
planetary perturbations are a plausible stirring mechanism 
if there are planets in the 1-30AU region, since the proximity 
of such planets to the disk means that their perturbations 
could affect the disk within the stellar age. For example, for 
a planet with an eccentricity of 0.1 to have stirred the disk 
at 30 AU within 4.6Gyr requires its mass to be > 150ap^^ (in 
Me), where api is its semimajor axis in AU (see eq. 15 of 
Mustill & Wyatt 2009). Such planets are not ruled out by 
the current radial velocity observations (see Fig. 9). 

4.6 Summary 

To conclude, we know that there are planets within lAU, 
and that there are planetesimals producing dust beyond 
30 AU, but the only thing we can say for certain about the 
1-30AU region is that there is a relative dearth of planetesi- 
mals within 30AU compared with the more distant regions. 
The presence of planets in the 1-30AU region would natu- 
rally explain the absence of planetesimals, but the emptiness 
of this region could also be explained by collisional erosion; 
i.e., it could have been full of km-sized planetesimals that 
were subsequently depleted through mutual collisions that 
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turned them into dust that was then removed by radiation 
forces. The latter explanation must also invoke a mechanism 
that stirs the disk (so that collisions are erosive). This would 
be readily explained if there were planets in the 1-30AU 
region, but could be inherent in the debris disk formation 
process. 



5 LOW-MASS PLANETS AND DEBRIS: IS 
THERE A CORRELATION? 

The coexistence of planets and debris around the 8th nearest 
G star prompted us to reconsider the lack of correlation 
between the two populations. In particular, we considered if 
earlier studies found no correlation (e.g., Greaves et al. 2004; 
Moro-Martm et al. 2007; Bryden et al. 2009) because at that 
time radial velocity studies were only sensitive to relatively 
massive planets. Here we suggest that a correlation does 
indeed exist for low-mass planets. 

5.1 Nearest 60 G star sample 

To assess the existence of a correlation with low-mass planets 
we considered a restricted but unbiased sample comprised of 
the nearest 60 G stars (see Phillips et al. 2010). The sample 
size was restricted to allow a thorough analysis of individual 
stars; an analysis of a larger sample is deferred to Moro- 
Martm et al. (in prep.). 

Out of this sample, 11 are reported in the literature to 
have planets, all of which were found in radial velocity sur- 
veys. We further subdivided this sample into those that have 
at least one massive planet, which we defined as a system 
with a planet more massive than > 95M^ (i-e., Saturn-mass) 
(hereafter called the high-mass planet sample), and those 
with only low-mass planets (hereafter called the low-mass 
planet sample). The high-mass planet sample comprises 5 
systems (HD146513, 47 UMa, /i Ara, 51 Peg, HD190360) 
and the low-mass planet sample 6 systems (HD20794, 61 
Vir, HD102365, HD69830, HD136352, HD38858) that are 
described further in Table 2. 

It may be questioned whether this division is appro- 
priate, given that longer observations may reveal systems 
currently classified as having only low-mass planets to have 
more massive planets at larger distances. However, the ra- 
dial velocities for 61 Vir show that there can be no planets 
more massive than Saturn orbiting within 6AU in this sys- 
tem (§4.1). So if more massive planets do exist they would 
have to be at larger separation, and the system would not 
have been classified as a high-mass planet system in this 
study. In other words, our classification of a low-mass planet 
system requires that there are planets detected above the 
blue line on Fig. 9, but that none of these is more massive 
than Saturn (or this would have been classified as a high- 
mass planet system), noting that the exact threshold will 
vary slightly between stars. 

Of the high-mass planet sample, all have been observed 
with Spitzer at 70 /im, and all but HD 190360 were also ob- 
served with Spitzer at 24 /xm (although fairly stringent lim- 
its do exist for this source from IRAS 25 /xm observations). 
None were found to have excess emission from circumstellar 
debris (Bryden et al. 2009). A detection rate of 0/5 is consis- 
tent with the number of detections that would be expected 



if the presence of debris is uncorrelated with the presence of 
planets, since nearby stars have a 14 — 16% detection rate 
with Spitzer (Trilling et al. 2008; Bryden et al. 2009). 

Note that if the Sun had been included in this sample 
Jupiter would fall above the detection threshold of radial 
velocity surveys, and so the Sun would have been classified 
as a high-mass planet system. Although we know the Sun 
has a debris disk, this is much lower in brightness compared 
with the disks known around nearby G stars (e.g.. Greaves 
& Wyatt 2010), and its thermal emission would not have 
been detectable in any current debris disk survey. This il- 
lustrates that the lack of debris disk detections around high- 
mass planet systems does not necessarily mean that these 
stars are completely clear of debris, since they could have 
disks below the detection threshold. Just as it was impor- 
tant to quantify what we mean by a system with high- and 
low-mass planets (see above, and Fig. 9), so it is important 
to also quantify what we mean by a system without debris; 
this is discussed further in §5.2.6. 



5.2 Debris around the low-mass planet sample 

In contrast to the high-mass planet sample, 4/6 of the low- 
mass planet sample are found to have debris. The debris 
toward three of these (61 Vir, HD69830, HD38858) was al- 
ready known about from Spitzer (Bryden et al. 2006; Beich- 
man et al. 2006). This is already a high detection fraction, 
because Poisson statistics show that the detection of > 3 
in a sample of 6 when the expected detection rate is 15% 
is a 6.5% probability event. This was not noted previously, 
because the planetary systems were only discovered since 
the Spitzer debris detections were reported. Excess emis- 
sion toward HD20794 was not reported from early analysis 
of the Spitzer observations (Beichman et al. 2006), but was 
detected by Herschel as part of the DEBRIS program, and 
this discovery is reported here for the first time. Below we 
provide more details on the debris toward members of the 
low-mass planet sample. 



5.2.1 HD20794 

At 6.0pc the G8V star HD20794 is the 5th nearest G star to 
the Sun; its age deduced from chromospheric Ca II H&K is 
6.2Gyr with 50% uncertainty (Vican 2012). A non-detection 
with Spitzer was reported in Beichman et al. (2006), but a 
reanalysis of Spitzer data, and a remodelling of the stel- 
lar photosphere, shows that at 70 /im the observed fiux 
of 106.7 d= 5.8 mJy is in fact significantly higher than the 
photospheric level of 80.9 mJy. The PACS observations also 
show a marginal excess at 70 /im of 97 =b 5.5mJy (the pho- 
tosphere is 84.2mJy in this bandpass), and a significant 
excess at 100 /im of 59.0 =b 4.8 mJy (compared with just 
41.0 mJy of photospheric fiux), both measured with a 15 arc- 
sec aperture. Although only marginal levels of excess were 
measured with MIPS at 24 /im (758.8 =b 7.6 mJy compared 
with a 741.5 mJy photosphere) and with PACS at 160 /im 
(20.6 =b 2.8 mJy from a PSF fit to the image, compared 
with a 15.8 mJy photosphere), these observations provide 
additional constraints that help pin down the temperature 
of the excess to 110 =b 35 K with a fractional luminosity of 
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Table 2. Nearest 60 G stars with planets 



Star name Planets (a, Mpi) 



Debris (r, /) 



References 



High-mass planet systems 

HD147513 (1.32AU, 400Me) 

(2.1AU, 840Me), (3.6AU, 180Me), (11.8AU, 550Me) 
(0.09AU, IIM©), (0.92AU, 170M©), (1.5AU, 560Me), (5.2AU, 600Me) 
(0.05AU, 160Me) 
(0.13AU, 18Me), (3.9AU, 500Me) 



47 UMa 
fi Ara 
51 Peg 
HD190360 



None detected 
None detected 
None detected 
None detected 
None detected 



1, 11 

2, 11 

3, 11 

4, 11 

5, 11 



Low-mass planet systems 

HD20794 (0.12AU, 2.8Me), (0.20AU, 2.5Me), (0.35AU, S.OMe) 

61 Vir (0.05AU, 5.3Me), (0.22AU, 19Me), (0.48AU, 24Me) 

HD102365 (0.46AU, 17Me) 

HD69830 (0.08AU, llM®), (0.19AU, 13Me), (0.63AU, 19Me) 

HD136352 (0.09AU, 5.5Me), (0.17AU, 12Me), (0.41AU, lOMe) 

HD38858 (l.OAU, 32Me) 



(19AU, 4 X 10-^) 
(29-350AU, 3 x 10"^) 
None detected"^ 
(lAU, 2 X 10-4) 
None detected"^ 
(102AU, 8 X 10-^) 



6, 12, 13 

7, 13 

8, 12, 13 

9, 14 

10, 12 
10, 12, 15 



"^^See Fig. 11 for constraints on disk luminosity, l=Mayor et al. (2004), 2=Gregory & Fischer (2010), 3=Santos et al. (2004), 4=Mayor 

& Queloz (1995), 5=Vogt et al. (2005), 6=Pepe et al. (2011), 7=Vogt et al. (2010), 8=Tinney et al. (2011), 9=Lovis et al. (2006), 
10=Mayor et al. (2011), ll=Bryden et al. (2009), 12=Beichman et al. (2006), 13=This work, 14=Beichman et al. (2011), 15=Lawler et 

al. (2009). 
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Figure 10. Spectral energy distributions for low-mass planet systems with debris: (left) HD20794, (middle) HD69830, (right) HD38858. 
Total fluxes are black points, open circles are predicted photospheric fluxes, grey points are photosphere subtracted fluxes. The photo- 
spheric model is shown with a blue line, and is a flt to data at wavelengths short enough that no excess is present. Longer wavelength 
fluxes are fltted as a modifled black body at the temperatures given in the text which is shown with a red line. The total flux from the 
models (photosphere + disk) is shown with a black solid line. 



4.3 X 10"^ (see Fig. 10)'^. The stellar luminosity of O.GGLq 
puts the dust at 5.2AU if it emits like a black body, but we 
estimate its location using the same /t = 1.9 factor derived 
for 61 Vir to predict its location to be around 19AU, not- 
ing that the uncertainty in the temperature, as well as that 
inherent in converting this to a radial location, means this 
location cannot be considered well constrained until the disk 
has been resolved.^ 



^ Note that the marginal excesses are plotted as upper limits on 
Fig. 10, but explain why the model flt for the excess flux on that 
flgure goes slightly under the two signiflcant detections. 
^ The 70 and 100 fim PACS observations show tentative evidence 
for extension at around the predicted scale, which is why fluxes 
were derived using an aperture at these wavelengths. 



5.2.2 HD 102365 

At 9.2pc the G2V star HD102365 is the 12th nearest G 
star to the Sun; its age deduced from chromospheric Ca 
II H&K is 5.7Gyr (Vican 2012). The previous non-detection 
with Spitzer (Beichman et al. 2006) is confirmed in our re- 
analysis of the data at 70 /xm (47.5 ± 4.5 mJy observed for 
40.1 mJy photosphere) and 24 /xm (363.1 ±3.6 mJy observed 
for 366.9 photosphere). This is further strengthened by a 
non-detection of excess in PACS images at both 100 /xm 
(15.9 ± 2.0 mJy observed for 20.3 mJy photosphere) and 
160 /im (3.0 =b 4.4 mJy observed for 7.8 mJy photosphere), 
where both measurements were made by fitting template 
PSFs to the images. This star has a binary companion at 23" 
projected separation, corresponding to an expected semima- 
jor axis of 297AU (Allen, Poveda & Herrera 2000). We fit the 
photometry available for the companion to estimate an ef- 
fective temperature of 2699K and note that its 0.7 mJy pho- 
tosphere is not detected in the 100 /im maps ( — 1.7 it 2.0 mJy 
was measured at the binary location). 
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5.2.3 HD69830 

At 12.6pc the G8V star HD69830 is the 22nd nearest G star 
to the Sun; its age deduced from chromospheric Ca II H&K 
is 5.7Gyr (Vican 2012). The excess toward this star has been 
analysed in great detail, since its infrared spectrum contains 
numerous spectral features that provide valuable informa- 
tion on the composition of the dust. The reader is referred 
to Beichman et al. (2011) for more details of the observa- 
tions, and here we summarise that the dust temperature 
puts it at ~ lAU, a location consistent with high resolu- 
tion mid-IR imaging and interferometry (Smith, Wyatt & 
Haniff 2009), and that its composition is similar to that of 
main-belt C-type asteroids in the Solar System. This loca- 
tion is at odds with the age of the star, since steady state 
collisional erosion of a planetesimal belt at lAU should have 
reduced the mass of such a belt to levels far below that ob- 
served (Wyatt et al. 2007a). This leads to the conclusion 
that the dust must be either transient (e.g., the product 
of a recent collision), or that it is replenished from a more 
distant planetesimal belt where longer collisional timescales 
mean that significant mass can survive (Wyatt et al. 2007a; 
Heng 2011). Far-IR observations have so far failed to find 
evidence for any outer belt (Beichman et al. 2011). 

5.24 HD136352 

The G2V star HD136352 is 14.8pc from the Sun; its age 
deduced from chromospheric Ca II H&K is 6.4Gyr (L. Vican, 
personal communication). No excess was reported toward 
this star with Spitzer at 70 /im (17.5 =b 5.1 mJy observed 
for 19.9 mJy photosphere; Beichman et al. 2006). This star 
was not observed by Herschel as part of the DEBRIS survey 
because of the high background level (Phillips et al. 2010). 

5.2.5 HD38858 

The G2V star HD38858 is 15.2pc from the Sun; its age de- 
duced from chromospheric Ca II H&K is 6.2Gyr (L. Vican, 
personal communication). Excess at 70 /im was detected 
with Spitzer (153.7 it 9.8 mJy observed for 15.0 mJy photo- 
sphere; Beichman et al. 2006), and later excess in the mid-IR 
spectrum was measured allowing the temperature to be de- 
termined and a dust location estimated at 12-29AU (Lawler 
et al. 2009). Like HD 136352, this star was not observed by 
Herschel as part of the DEBRIS survey because of the high 
background level (Phillips et al. 2010). For consistency with 
our interpretation of other sources, we repeat the analysis of 
its excess spectrum (see Fig. 10), and use the inferred dust 
temperature of 50K to estimate a location from this O.83L0 
star as 102AU (again using the /t = 1.9 factor inferred for 
61 Vir). This size is consistent with the marginal extension 
of the MIPS 70 /xm emission that can be fitted by a ring of 
size around 100 AU size (section 5.2.3 of Lawler et al. 2009). 

5.2.6 Summary 

The detection rate of 4/6 is much higher than that seen in 
the background population. Although the analysis included 
observations with Herschel, in fact all of the detections were 
present in the Spitzer data, and so it may be reasonable to 
compare this detection rate with the 15% of Sun-like stars 



found to have excess by Spitzer. Poisson statistics shows that 
detecting > 4 from a sample of 6 for an expected detection 
rate of 15% is a 1.4% probability event. If instead the 4/6 
rate was compared with the detection rate of disks found 
by Herschel in an unbiased population of Sun- like stars^, 
which is around 20-25%, a rate this high would have been 
up to a 6% probability event (for a 25% detection rate). 
Thus, although this analysis is necessarily subject to small 
number statistics, it appears to be very unlikely that the 
low-mass planet sample has as low a detection fraction as 
other Sun-like stars. As one of the two low-mass planet sys- 
tems without detected debris is a wide binary system, the 
high disk fraction would be further enhanced if the sam- 
ple was restricted to single low-mass planet systems, though 
having a wide binary companion does not necessarily pre- 
clude the presence of debris (Trilling et al. 2007; Rodriguez 
& Zuckerman 2012). 

Note that we are only claiming that the debris incidence 
rate is higher around stars with low-mass planets (that are 
detectable in current radial velocity surveys), not that all 
low-mass planet systems have detectable debris. Indeed this 
appears unlikely given that low-mass planetary systems may 
be present around 30-50% of stars (Mayor et al. 2011), and 
that only 15% of stars have debris that was detectable with 
Spitzer. Also we cannot comment on any correlation with 
planets closer to Earth-mass that would be below the detec- 
tion threshold of radial velocity surveys (e.g.. Fig. 9), though 
such planets may be common (e.g., Melis et al. 2010; Borucki 
et al. 2011). 

Likewise it is important to clarify that a detectable de- 
bris disk in this context implies that the fractional luminos- 
ity and disk temperature lie above the appropriate lines on 
Fig. 11. Although the detection thresholds are similar for 
the different stars in the low-mass planet sample, there is a 
dispersion in sensitivity to large cold disks (< lOOK), and it 
should be acknowledged that a disk with the same fractional 
luminosity and temperature as that of HD 20794 would not 
have been detectable around all of the sample. Neverthe- 
less, the same would be true of disks that could have been 
detected around the rest of the nearest 60 G star sample. 

We defer a detailed consideration of the planet correla- 
tion to a later paper (Moro-Martm et al., in prep.), but note 
that the discovery of a debris disk around an M star with 
low-mass planets could suggest that this trend also continues 
to lower mass stars (Lestrade et al., in prep.). 

5.3 Implications for planet formation 

The correlation summarised in §5.2.6 suggests that, even 
if not directly connected, the formation of a system that 
was conducive to the formation of low-mass planets was 
also conducive to large quantities of planetesimals remain- 
ing in the system after several Gyr of evolution. This may 
be connected to a different correlation, in which stars with 
high-mass planets tend to be metal-rich whereas those with 
low-mass planets have no bias in metallicity (Sousa et al. 

^ For example, 38 of the nearest 60 G star sample have been 
observed by Herschel as part of either the DEBRIS or DUNES 
key programmes (Matthews et al. 2010; Eiroa et al. 2010), of 
which 8 have disk detections. 
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Figure 11. Parameter space for which excess emission from a 
debris disk would have been discovered in existing observations 
of the stars in the low-mass planet sample. The limits are de- 
rived by assuming that the star has a disk with a given frac- 
tional luminosity (Lir/L^) and temperature (which is translated 
into disk radius using /t = 1.9), and determining whether that 
emission would have been detectable in the observations at any 
wavelength at which the star was observed using eq. 8 and 11 
of Wyatt (2008) with the modification of Lestrade et al. (2009) 
for disks extended beyond the beam. The limits at different disk 
temperatures are dominated by different wavelengths (e.g., tem- 
peratures above 170K are limited by MIPS 24 /im observations, 
with longer wavelengths constraining lower temperatures). The 
limits for different stars in the sample are plotted with different 
colours, and the black dots denote the parameters derived for the 
disks detected around four of the sample. 



2011), and neither do stars that host debris disks (Greaves 
et al. 2006; Maldonado et al. 2012). 

Ideally the debris disks in low-mass planetary systems 
would share some characteristics that could both point to 
a reason for the correlation between these phenomena, and 
help identify systems as having low-mass planets through 
study of their debris disks. However, Table 2 shows that de- 
bris disks in low-mass planetary systems have a range of radii 
( 1-350 AU) that spans the range of known debris disks. Thus, 
planetesimals in low-mass planetary systems could either be 
at large distances from the planets (like 61 Vir, HD20794 
and HD38858), or some of those planetesimals could be in 
close proximity to the planets, at least at some point in their 
evolution (like HD69830). 

Why might such a correlation exist? Here we consider 
the two scenarios that have been proposed for the formation 
of low-mass planet systems: (i) that the planets formed fur- 
ther out but then migrated in, perhaps gaining mass as they 
migrate (Alibert et al. 2006, Terquem & Papaloizou 2007, 
Ida & Lin 2010), or (ii) that they formed in situ in which 
case significant solid mass (50 — lOOM®) would have to be 
concentrated within lAU (Hansen & Murray 2011). 

Although not without issues in comparison with exo- 
planet statistics (McNeil & Nelson 2010), formation scenario 
(i) seems to be perfectly suited to leaving large quantities of 
outer debris. Payne et al. (2009) showed that the migration 
of the planets through a planetesimal disk is an inefficient 
process, in the sense that not all of the planetesimals are 



accreted or ejected by the migrating planets. Rather, many 
planetesimals are left stranded in the region through which 
the planets recently migrated. As long as there were no more 
distant planets that subsequently migrated in thus depleting 
this region, these planetesimals could be dynamically sta- 
ble over long timescales, but would predominantly deplete 
through mutual collisions. 

Formation scenario (i) could also be applied to the for- 
mation of close-in massive planets. However, in this case the 
stronger gravitational perturbations of the migrating planets 
would result in the ejection or accretion of a larger fraction 
of planetesimals (though not necessarily all of them, e.g., 
Raymond et al. 2006). Moreover the migration of multiple 
high-mass planets may result in a dynamically unstable con- 
figuration, unless this can be prevented by resonance lock- 
ing, in which case the ensuing instability could deplete any 
remaining planetesimals (e.g.. Gomes et al. 2005; Raymond 
et al. 2011). 

Formation scenario (ii) may also be conducive to form- 
ing an outer planetesimal belt, depending on how efficiently 
the solid mass is accreted onto planets. If even a small frac- 
tion of the planetesimal mass is scattered onto highly ec- 
centric orbits, it would form a mini-scattered disk of plan- 
etesimals with pericentres in the vicinity of the planets but 
apocentres extending out to large distances. If this happened 
during the protoplanetary disk phase, interaction with solid 
or gaseous material in the outer planetary system could then 
pull the pericentres of scattered planetesimals away from the 
planets (i.e., through gas drag or dynamical friction). This 
could prevent the planetesimals from being ejected and leave 
them in a long-lived extended disk configuration. 

Again, application of scenario (ii) to the formation of 
massive planets would be less favourable to the formation of 
massive outer debris disks. This is because planetesimals are 
more readily ejected by more massive planets (e.g., Tremaine 
1993), on timescales that would be shorter compared with 
the timescale for the outer disk to decouple the planetesimals 
from the planets' gravitational influence. 

It thus seems plausible that high levels of outer debris 
could be a natural outcome of systems that form only low- 
mass planets. 



6 CONCLUSIONS 

Observations from the DEBRIS Key Programme are pre- 
sented that resolve the structure of the debris disk around 
the exoplanet host star 61 Vir (§2). Modelling shows that 
the dust extends from 30 AU out to at least lOOAU in a 
nearly edge-on configuration (§3). There is likely little inter- 
action between the disk and the known planets, which are 
at < 0.5AU. The lack of planetesimals in the < 30 AU region 
could be explained by the existence of planets in this region, 
but the depletion can also be explained by collisional erosion 
(§4). Considering a sample of the nearest 60 G stars there is 
an emerging trend that stars which, like 61 Vir, only have 
low-mass planets, are more likely to have detectable debris 
(§5). We attribute this trend to the fact that the formation 
processes that make low-mass planets are likely to also result 
in large quantities of distant debris. 
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